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ABSTRACT 

As a planet transits the face of a star, it accelerates along the line-of-sight. The changing delay in 
the propagation of photons produces an apparent deceleration of the planet across the sky throughout 
the transit. This persistent transverse deceleration breaks the time- reversal symmetry in the transit 
lightcurve of a spherical planet in a circular orbit around a perfectly symmetric star. For "hot Jupiter" 
systems, ingress advances at a higher rate than egress by a fraction ~ 10~''-10~'^. Forthcoming space 
telescopes such as Kepler ot COROT wiW reach the sensitivity required to detect this asymmetry. The 
scaling of the fractional asymmetry with stellar mass and planetary orbital radius a as oc Mi, /a?' 
is different from that of the orbital period as oc (M^,/a^)~^/^. Therefore, this effect constitutes a 
new method for a purely dynamical determination of the mass of the star, which is currently inferred 
indirectly with theoretical uncertainties based on spectral modeling. Radial velocity data for the reflex 
motion of the star can then be used to determine the planet's mass. Although orbital eccentricity 
could introduce a larger asymmetry than the light propagation delay, the eccentricity is expected to 
decay by tidal dissipation to negligible values for a close-in planet with no perturbing third body. 
Future detection of the eclipse of a planet's emission by its star could be used to measure the light 
propagation delay across the orbital diameter, 46.7(a/7 x 10^^ cm) seconds, and also determine the 
stellar mass from the orbital period. 

Subject headings: planetary systems, techniques; photometric 



1. INTRODUCTION 

The population of known extrasolar planets which 
transit the face of their parent stars has been grow- 
ing steadily in recent years. It currently includes 
HD 209458b lICha rbonnea ^ireraLl I2OO0I) OGLE-TR- 
56b llTorre s et al."2004||). 0GLE-TR-113b fKonack i et all 
l20n4aD. GLE-TR-132b (||loutou et al. 2004). TrES-1 
ijAlonso et al. 2004), O GLE-TR-lll b (Pont ct al. 2007) 
and OGLE-TRlOb (,Konacki et all 12004^ . For HD 
209458b , the transit lightcurve was observed with the 
Hubble Space Telescope (HST) to be symmetric around 
its centroid to an exquisite photometric precision of 
~ 10~^ magnitude per data point for a few hundred 
data points tracing both ingress and egress ijBrown et al.l 
l2f)0l . 

In this Letter we show that transit lightcurves must 
have a time-reversal asymmetry even if the star and the 
planet are perfectly symmetric and the orbit is circular. 
The asymmetry originates from the persistent accelera- 
tion of the planet towards the star during the transit. For 
a circular orbit, the planet moves towards the observer at 
the beginning of the transit (ingress) and away from the 
observer at its end (egress). This net acceleration intro- 
duces a change in the relative rate by which ingress and 
egress advance in the observer's frame of reference. The 
effect simply follows from the unsteady change in the 
propagation delay of photons which distorts the trans- 
formation of time between the planet and the observer. 
The variation rate of the delay is changing most rapidly 
at the middle of the transit when the full Newtonian ac- 
celeration vector points straight along the line-of-sight. 
In the following sections, we derive the magnitude of the 
resulting lightcurve asymmetry (§2) and discuss its im- 
plications (§3). 



2. APPARENT TRANSVERSE DECELERATION DUE TO 
PROPAGATION DELAY OF PHOTONS 

The observed arrival time of photons tobs is given by 
the time they left the planet t plus the propagation delay 
over the distance of the planet D, 

tobs = t + D/c. (1) 

The ^-derivative of equation Q gives 

dtobs = rfi (1 + ^) , (2) 

where v\\{t) — dD/dt is the velocity of the planet along 
the line-of-sight. Since t;|| changes from negative to posi- 
tive during the transit of a planet in a circular orbit, it is 
now obvious that ingress would advance at a faster rate 
(i.e. a shorter d^obs per dt interval) than egress. Figure 
1 illustrates schematically this generic behaviour. Note 
that equation |2Jl is accurate to all orders in |v/c| but 
we will keep only first order terms in subsequent deriva- 
tions. Even though the planet is not emitting the ob- 
served photons, the timing of its occultation is dictated 
by the instant at which stellar photons graze its outer 
surface; this timing is distorted by the same propagation 
delay effect that would exist if the grazing photons were 
emitted by the planet itself rather than the star since 
the propagation history of the photons before they graze 
the planet's boundary is irrelevant. The transverse spa- 
tial coordinates of the planet orbit, xj^, are the same in 
the observer and planet reference frames. The appar- 
ent transverse velocity of the planet is therefore differ- 
ent from it s actual transverse velocity, vj^ ~ (dx^/rfi) 
ijRvbicki fc Lightman 1979jl , 

•° ~ dtobs \dtobsJ dt \ c J ' 
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Fig. 1. — Schematic illustration of the propagation delay effect. 
The transiting planet is moving towards the observer (d|| < 0) dur- 
ing ingress and away from the observer (d|| > 0) during egress. As 
a result, the observed time interval of partial eclipse is shorter at 
ingress relative to egress (see Eq. 1^. This breaks the time re- 
versal symmetry of the observed lightcurve and introduces a frac- 
tional difference of <Sd||/c in the temporal slope of the lightcurve 
between ingress and egress, where (5i)|| = iiyr is the net gain in 
the planet's line-of-sight velocity over the transit duration t. For 
close-in planets, the fractional asymmetry in the slopes is of order 
S ~ 10-*-10-3 (see Eq. |7}. 



where the transverse position vector x.± corresponds to 
angular coordinates on the sky times the distance to the 
planetary system. Throughout our discussion, the terms 
parallel (||) or transverse (_L) are relative to the line-of- 
sight axis that starts at the observer and goes through 
the center of the stellar image. We focus on the time 
when the occulting planet crosses this spatial point of 
symmetry. 

Taking the tobs~derivative of both sides of equation (jJJ 
and keeping terms to leading-order in |v/c|, we get two 
new contributions to the difference between the observed 
and Keplerian values of the transv erse acceleration of the 
planet across the sky l|Loebll2003|) . 
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(4) 



where Vobs = i'P^/dtl^^s) is the observed acceleration 
and V = {(Px/dt'^) is the actual Keplerian acceleration 
of the planet. The last term on the right-hand-side of 
equation Q implies that the apparent transverse accel- 
eration of the planet gets a contribution from its Keple- 
rian acceleration along the line-of-sight, This term is 
the source of our effect. 

For simplicity, we assume that the orbital plane is 
viewed edge-on by the observer; this geometry is exact 
for transits that cross the center of the star ("central 
transits" ) and is a very good approximation more gener- 
ally as long as the orbital radius is much larger than the 
radius of the star. For a circular planetary orbit, Newto- 
nian dynamics implies no transverse acceleration of the 
planet at the spatial center of symmetry of the transit 
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(5) 



Here Af* is the stellar mass and a is the planet's orbital 
radius. However, equation Q implies that at the same 
time there would be an apparent transverse deceleration 
of the planet in the observer's frame of reference 



c " 



(6) 



which is directed opposite to v^. This persistent ap- 
parent deceleration implies that the planet would cross 
ingress and egress at different apparent rates. Thus, even 
if the stellar image is perfectly symmetric, the orbit is 
circular and the planet is perfectly spherical, there is an 
inherent asymmetry in the transit lightcurve due to the 
nearly steady value of = GMi^ja? during the transit. 
Newtonian dynamics alone predicts a non-vanishing vj^ 
as soon as the planet moves away from the transit center, 
but for a circular orbit this deviation would maintain the 
time-reversal symmetry between ingress and egress. 

The fractional asymmetry in the transit lightcurve is 
of order the fractional change in v±^ohs over the transit 
duration, since the rates by which ingress and egress pro- 
ceed are proportional to w^^obs- For a total transit dura- 
tion T and an orbital period T = 2TTa/v±^ 3> t, the slope 
of the initial drop and final rise in the transit lightcurve 
will differ by a fractional amplitude 
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For a "hot Jupiter" like HD 209458b , the ingress phase 
< 0) would proceed at a rate that is higher by 
fractional amplitude of ~ 10"'* than the egress phase 
> 0). Planets that are closer in by a factor of a few 
could produce an asymmetry of up to (5 ^ lO^'^; OGLE- 
TR- 56b for which a = 3 5x 10" cm provides S = 3x IQ-'^ 
[see IGaudi et al.l l)2005j) for a compilation of all known 
transit systems and the prospects for detecting others]. 

The observed time tohs is a slightly distorted version of 
the time axis t along which the lightcurve is symmetric. 
We may write tobs = where e{t) — {t/2T)S and we 

shifted i = to be at the transit centroid. The observed 
photon flux, F{tohs), corresponds to the flux at an undis- 
torted time t « iobs[l while F{—tohs) corresponds 

to i w — iobs[l — ^i^t)] — ~^obs[l + For small devi- 
ations, we may expand the photon flux as a function of 
time to leading order. Fit + At) « F{t) + {dF/dt)\tAt. 
Since e{t) = -e{-t), F{t) = F{-t), and {dF/dt)\t = 
— {dF/dt)\-t, the time-dependent photometric asymme- 
try of the lightcurve is given by 



... X f ftobs) - -F(-tobs) 

F{^) 

dhiF 



(8) 
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where we have set tobs = at the transit centroid (for 
which dF /dt ohs = 0) and kept terms to leading order in 



3 




oba 



(hours) 



Fig. 2. — Predicted fractional deviation from time-reversal 
symmetry, Ap(tchs)j for t he observed transit lightcurve of HD 
209458b (Brc moret al.,:2QQll) . based on Eqs. Q and (8). The verti- 
cal axis is proportional to (M*/l.lMQ)(a/7x lO^^cm)"^. The four 
breaks in the curve result from the four points of contact between 
the projected planet boundary and the limb of the star. 



|e| <C 1. Given a preliminary transit lightcurve, -F(iobs)i 
and an approximate knowledge of the system param- 
eters, it is possible to predict the lightcurve asymme- 
try under the assumption of a circular orbit. Figure 2 
shows the expected A p(^nbs) for HP 2094 58b based on 
the lightcurve data in Bro wn et al.l l)2001|) . For a tran- 
sit depth of 2 X 10^^ magnitudes produced by a close-in 
planet, the net photometric asymmetry would typically 
amount to a few percent of 6 or equivalently to a peak 
value of |Ai?| in the range of (O.l-l)xlO^^. The photo- 
metric sensitivity of 10~* magnitudes per data point 
achieved with HST for HD 209458b ( Brown ct al. 2001) 
provides a net sensitivity to a time-reversal asymmetry 
of |Ai.| - 10-5(iV/102)-i/2^ where N ~ 10^ is the total 
number of independent data points during ingress and 
egress. Our effect appears to be on the borderline of be- 
ing detectable with existing techniques. A statistically 
significant (> Scr) detection of the asymmetry in equa- 
tion lO for HD 209458b would require an additional or- 
der of magnitude improvement in sensitivity, which is 
achievable with future space missions such as Kepler^ 
or COROT ^ and only marginally with the existing in- 
strument MOST A simple way to reach the required 
sensitivity would be to maintain the existing photomet- 
ric precision per data point and increase the number of 
data points from ^ 10^ up to ~ 10'* (e.g. by increasing 
the number of HD 209458b transits from the 4 observed 
by HST to ~ 400) over a period of several years. 

The orbital period scales as oc (M^/a^)~*/^ while the 
fractional asymmetry in equation Q scales as oc Mi^ja? . 
Thus, detection of the lightcurve asymmetry would al- 
low to determine M^, and a separately. Additional data 
on the reflex motion of the star would then provide the 

t http:/ /www. kepler.arc.nasa.gov/summary.html 
^ http:/ /serweb. oamp.fr/projets/corot/ 
3 http://www.astro.ubc.ca/MOST/ 



planet mass, Afp, based on dynamical information alone. 
Currently, the stellar mass, M*, is inferred indirectly 
based on theoretical fitting of the stellar spectrum which 
is subject to modeling unce r tainti es at the level of > 10% 
[see Fig. 3 in I Alonso et"al] ()2n04D ]. 

If the planetary orbit is nearly circular but has a finite 
eccentricity, e ^ 1, then it is easy to show that at the 
transit center 

w^^-esm^- — (9) 

to leading order in e, where is the angle between 
the shortest radius vector of the orbital ellipse (along 
the apsidal line) and the line to the observer. As long 
as |e sin 5*1 < |wx/c| < 10"'^, the eccentricity effect is 
smaller than the propagation delay effect in equation (|SJ). 
In the absence of perturbers, any initial orbital eccen- 
tricity is expected to decay exponent ially through tidal 
dissip ation on an e-folding timescale l|Goldreich fc Soteil 
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where the quality parameter Qp (~ 10^ for Jupiter) is in- 
versely proportional to the dissipation rat e in the p lanet's 
interior (loannou & Lindzcn 1993; Qgilvi e^ Linl l2004'). 
and Rp is the planet's radius. Aside from the unknown 
Qp, the parameter values in equations Q and (|10|l were 
chosen to match HD 209458b . Since the orbital cir- 
cularization timescale is shorter by up to two orders 
of magnitude than the age o f planetary sy stems such 
as HD 209458 5 Gyr, see iMazeh et all l|2000j) and 
.Cody fc Sas selov (2002^], the eccentricity-driven asym- 
metry is expected to diminish for a close-in planet unless 
a third body pumps its orbital eccentricity. Incidentally, 
such a perturber was hypothe sized as the driver of the 
inflated radius of HD 209458b l)Bodenheimer et alJl2003() 
but its existence has not been demonstrated yet. 

For close-in planets such as HD 209458b , our effect 
is larger than a different source of lightcurve asymme- 
try that was already di scussed in the literatu re, namely 
the planet's obliquity l|Hui & Seageil l2002bH . This is 
because close-in planets are expected to rotate slowly 
and possess a small projected obliquity as a result of 
strong tidal locking with their orbital revolution. Their 
rotation axis is expected to be normal to their orbital 
plane which is viewed nearly edge-on [in the case of HD 
209458b , the orbita l inclination angle is i = 86.° 1 ± l.°6 
1 1200(1 ] 



fMazeh et alJ 120001) ]. Another source of asymmetry is 
rotation of the star. Ty pical projected rotation speeds of 
Vr„t s in i 4 km s"* l|Mazeh et all 12001 iQueloz et alJ 
I 20QQ |) would produce a Doppler offset fLoeb & Gaudi' 
2003) between the flux emitted by the approaching and 
receding sides of the star of order ~ (wrot/c) = 1.3 x 10"^ 
[the pho tometric an alog of t he Rossiter- McLaug hlin ef- 
fect; see IQueloz e~a l. ( 2000'l: lOhta et al.l l|2004Jl ]. or an 
occultation contrast of ^ 10~^ magnitudes which is much 
smaller than the effect considered here. If the stellar ro- 
tation axis is significantly misaligned with the normal 
of the orbital plane and the transit is not central, then 
the rotation-induced oblateness of the stellar image could 
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generate a fractional ((5-equivalent) asymmetry of less 
than i(w^Qji?*/GA/^) ~ 5 x 10^^, which is again well be- 
low our effect for HD 209458b where the misalignment 
angle must be small (Queloz ct al. 2000). 

The propagation delay relation between tobs and t in 
equation can be easily incorporated into a computer 
program that searches for the best-fit Keplerian orbit 
under the constraints of a given data set. The delay- 
corrected Keplerian fit would provide new dynamical con- 
straints on the planetary system. Such a fit would involve 
the same number of free parameters as the standard Ke- 
plerian fit. 

3. DISCUSSION 

We have shown that the propagation delay of light in- 
troduces an apparent transverse deceleration of a planet 
on the sky during its transit across the face of its par- 
ent star (Eq. 0). This persistent deceleration breaks 
the time-reversal symmetry of the transit lightcurve for 
a spherical planet in a circular orbit around a spher- 
ical star. Throughout the transit, the planet velocity 
along the line-of-sight, t)|| , is changing at a nearly steady 
rate, GMi^/a^. This produces a steady change in the 
transformation of time to the observer's frame (see Eq. 
|5J). The net change in v\\ between ingress and egress, 
5v\\ = W||T, introduces a fractional difference of magni- 
tude 5 = 5v\\/c in the slopes of their lightcurves. It is 
possible to search for this difference in slopes by fold- 
ing the lightcurve over its centroid. For close-in planets, 
ingress should typically proceed at a rate that is faster 
by a fractional amplitude 5 ~ 10~'*-10~'^ than egress 
(see Eq. [Tj). This level of asymmetry will be detectable 
with forthcoming space telescopes such as COROT or 
Kepler (:,Borucki et a l. 2003), which are scheduled for 
launch within 2-3 years. Because this asymmetry has 
a unique scaling with stellar mass and orbital radius 
(oc Mi./a^), its detection together with the reflex motion 
of the star will allow to determine the star and planet 
masses as well as the orbital radius using purely dynam- 
ical data. This method evades the theoretical uncertain- 
ties inherent in the existing approach for determining 
stell ar masses based on mo deling of spectroscopic data 
[e.g. ICodv fc Sassekwl |2002)]. 

Orbital eccentricity could induce a stronger asymmetry 



in the lightcurve but is expected to decay exponentially 
to negligible levels through tidal dissipation for close-in 
planets like HD 209458b , unless it is being pumped by 
the gravitational perturbation of another planet. 
Similarl y to other transit timing residuals 



fMira,lda,-E^^ UHl lHolma,n fc Mun^ UHl 
lAgol et al.l 12004 ). our effect would be contaminated 



by noise from inhomogeneities on the face of the 
star and would compete against other small effects 
involving asymmetries from the oblateness of the planet 
ijHui fc Seager 2002a b) or the rotation of the star. 
Additional special-relativistic or general-relativistic 
effects are of the order of ~ {v/c)'^ or ~ ^/c^, or smaller, 
where is the gravitational potential produced by the 
star {(/) ~ u^); these corrections are orders of magnitude 
smaller than the propagation delay effect discussed here. 

Finally, we note that a change of the opposite sign in 
W|| occurs when the planet goes behind the star. In this 
case, ingress would be slower than egress. There is no net 
change in the orbital period over a full closed orbit. How- 
ever when the planet enters its own (secondary) eclipse 
by the star, the photons it emits will be delayed relative 
to primary eclipse (the transit) by the difference in emis- 
sion times plus the light travel time across the orbital di- 
ameter, as implied by equation For a circular orbit, 
the time interval between the centroids of the primary 
and secondary eclipses would be longer by ST1/2 = 2a /c 
than half of the full orbital period. Future detection of 
a planet's infrared emission would then allow to deter- 
mine the orbital radius from the light propagation delay 
of STi/2 = 46.7(a/7 x 10" cm) seconds. The fuh or- 
bital period T would then yield the stellar mass through 
M^, = a^/[G'(r/27r)^]. An eccentricity could change the 
lightcurve history of the illuminated planet and in par- 
ticular make the time from primary eclipse to secondary 
eclipse different from the time back to primary eclipse. 
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